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Introduction. Cylindrical and spherical shells are extensively used in engineering. They face internal and/or external 
pressure and heat. Stresses and strains distribution in elastoplastic shells has been studied by many scientists. Numerous 
works involve the use of the von Mises yield conditions, maximum shear stress, maximum reduced stress. These condi- 
tions do not include the dependence on the first invariant of the stress tensor and the sign of the third invariant of the 
stress deviator. In some cases, it is possible to obtain numerical-analytical solutions for stresses, displacements and de- 
formations for bodies with spherical and cylindrical symmetry under axisymmetric thermal and force action. 

Materials and Methods. The problem on the state of a thick-walled elastoplastic shell is solved within the framework of 
the theory of small deformations. A plasticity condition is proposed, which takes into account the dependence of the 
stress tensor on three independent invariants, and also considers the sign of the third invariant of the stress deviator and 
translational hardening of the material. A disconnected thermoelastoplastic problem is being solved. To estimate the 
stresses in the region of the elastic state of a spherical shell, an equivalent stress is introduced, which is similar to the 
selected plasticity function. The construction of the stress vector hodograph is used as a method for verification of the 
Stress state. 

Results. The problem has an analytical solution for linear plasticity functions. A solution is obtained when the strength- 
ening of the material is taken into account. Analytical and graphical relationships between the parameters of external 
action for the elastic or elastoplastic states of the sphere are determined. For a combined load, variants are possible 
when the plastic region is generated at the inner and outer boundaries of the sphere or between these boundaries. 
Discussion and Conclusions. The calculation results have shown that taking into account the plastic compressibility and 
the dependence of the plastic limit on temperature can have a significant impact on the stress and strain state of a hollow 
sphere. In this case, taking into account the first invariant of the stress tensor under the plasticity condition leads to the 
fact that not only the pressure drop between the outer and inner boundaries of the spherical shell, but the pressure values 
at these boundaries, can vary within a limited range. In this formulation of the problem, when there is only thermal 
action, the hollow sphere does not completely pass into the plastic state. The research results provide predicting the 
behavior of an object (a hollow sphere) that experiences centrally symmetric distributed power and thermal external 


influences. 
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Introduction. The solution to the problem of a thick-walled spherical shell experiencing different external 
influences is given in the monographs [1, 2] and a number of scientific papers on the theory of elasticity, plasticity, and 
thermoelastic plasticity [3-9]. Usually, the case is considered when the loading process is simple. 

The problem of a thick-walled spherical shell is one of the simplest elastoplastic problems when the fields of 
external actions and internal parameters have central symmetry. Due to the central symmetry in the plastic region, the 
regime of complete plasticity 1s performed. For an ideal plastic body, the problem is statically definable, which allows it 
to be solved under any plasticity conditions. In the monograph [1], the most complete solution to the spherical shell 
problem is given, when the dependence of the plastic limit on temperature is not taken into account, and the plasticity 
condition does not depend on the first invariant of the stress tensor and the sign of the third invariant of the stress 
deviator. The cases of only thermal and combined loading are considered, when the temperature on the walls of the 
sphere is set, the pressure on the inner walls is set, and there is no pressure on the outer wall. In [10-13], 
thermoelastoplastic state of various objects was studied, and in [3-9], the process of thermal loading and unloading of a 
sphere free from external forces and a hollow sphere for the Tresca condition was considered with account for the 
dependence of the plastic limit on temperature. The solution to this and similar problems is of interest since it is 
possible to obtain an analytical or partially analytical solution for various mathematical models. An analytical solution 
can be obtained through selecting piecewise linear plasticity functions [11, 13]. Mathematical modeling of objects 
enables to predict their state and behavior depending on the values of the initial parameters [14, 15]. 

Materials and Methods. Problem Statement. We consider the problem of a thick-walled spherical shell (a 


hollow sphere) experiencing centrally symmetric external influences: pressure p, on the outer wall at p=b5 and 
pressure p, on the inner wall at p=a. The thermal effect on the sphere is also considered: temperature 7 is 
maintained at the boundary p=a, temperature 7) is maintained at the boundary p=D. It is assumed that the sphere 


exhibits elastic and plastic properties. The desired state parameters at each point of the sphere are the components of the 
stress tensor, the components of the strain tensors, and the displacement vectors. In the elastic state region, the elastic 
deformations are complete (there are no residual deformations). 

Basic Ratios. All relations are reduced to a dimensionless form. The outer radius of the sphere b is selected as a 
length scale. All values having the stress dimension are assigned to the plastic limit under uniaxial tension k . The scale 
unit for temperature is | C. 


Due to the specified symmetry of external actions, in the spherical coordinate system p,0, of the matrix, the 


components of the stress and strain tensor will have the form: 


o, 0 0 é, O 0 
(s)=| 0 o, OF}, (e)=| 0 &s O 
0 0 Oo, 0 0 6, 


In this case, the equalities 6, =0,, €&, =, are fulfilled. 


o 9 
If the plasticity functions do not depend on the first invariant of the stress tensor and the sign of the third 
invariant of the stress deviator, then, when solving the sphere problem, the plasticity functions will be reduced to the 


form: 
f =|9,-9, Ek. (1) 


Consider the plasticity condition: 


e((G,-0e,) 4 200,06.) 2 


0.,0,, &,&) J= 
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When the parameters have the values: ¢=0, n=1,6=0,a =0, m=1, u=0,k =k,, the condition (2) implies the 
condition (1). In Fig. 1, the plane o,,0, shows the plasticity curves determined from the formula (2) for different 


values of the numerical coefficients in the plasticity function. 





a) b) 
Fig. 1. Plasticity curves: a) for parameters: ¢ =0,2; w= 2; 1 =0,5; 6 =0; m =3; k =1; solid line a =0,5; 


dotted line a@ = 0; 5) for parameters: o =0;¢=0;5=0;4 =1 


The results presented in Fig. 1 show that when the first invariant of the stress tensor is taken into account, the 
radial and circumferential stresses can vary in a limited range when the point of the sphere is in an elastic state. 
Accordingly, the pressure on the boundaries of the sphere should also be limited. When the first invariant in the 
plasticity condition is not taken into account, the elastic state is possible for any pressure value at the boundaries of the 
sphere, but the pressure drop is limited Ap = p, — p,. Taking into account the sign of the third invariant of the stress 


deviator, as noted above, affects the values of the plasticity limits. 
If the values of the state parameters 0,0, determine the point of the region bounded by the plasticity curve, it is 


assumed that the defining equations connecting stresses and deformations are the relations of the Duhamel-Neumann 
law [1, 2]: 
Ee, =(l-v)o,-vo,+ kal, Ee, =0,—-2vo,+E£a?, (3) 
where the Young's modulus F' and the Poisson's ratio v are constants. 
If the state parameters 6,,0, determine the points on the plasticity curve, then an additive representation of the 
total deformations in terms of reversible and irreversible deformations is assumed: 
E,=E, +, € =E +E. (4) 


Complete deformations are determined through displacements from the formulas: 





u du 
2) oo (5) 
p dp 
Complete deformations are bound by the condition of compatibility of deformations: 
d 
ig “8 4+¢,-8 =0. (6) 
r 


Increments of irreversible deformations are related to stresses by the normal law: 
ds’ 
de 2, (7) 
Of/0c, of /0o, 
The relation (7) is generally non-integrable when choosing nonlinear plasticity functions [16]. In the quasi-static 
approximation, the stresses must satisfy the equilibrium equation: 





do, 
fe Ip +2(6,—-0,)=0. (8) 


Equivalent Stress. The equivalent stress is the convex isotropic scalar functions of the stress tensor. In special 
cases, the term “equivalent stress” is synonymous with other terms, for example, “stress intensity” [17]. In this paper, 
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the equivalent stress coincides with the plasticity function. In this case, the equivalent stress will not have a 
discontinuity at the elastic-plastic boundary. 
Temperature Field. The temperature field in the sphere is found from the solution to the boundary value 
problem [1]: 
dT dT 
pt 0, 
dp dp (9) 
T la = ie T p= fie 





The solution to the problem (9) is presented in the form: 
Pe a re. (10) 
(b—a)\p 
Elastic Area. In the region of the elastic state of a hollow sphere, the formulas for stresses have the form: 
BAX B KX abEaAT 
0-Ar -—, 0, =A-_2-_, ha 
pp 2p 2p (l1—v)(b—a) 
Plastic Area. We select the conditions (1). Let us consider the case of only the thermal effect (10). Then, the 
plastic region will originate at the inner boundary of the shell under the condition [1]: 
Ea 2(a° +ab+b’)k 
p=" jarpp =A ter 
(l—v) b(a+2b) 


Denote by c, — the radius of the elastic-plastic boundary p =c,. During the loading process, when B > B,, the 











plastic area a<p<c, increases. When the condition (1) is selected, the stresses in the region a<p<c, are calculated 
from the formulas: 

o! =2KkIn(p/a), 6,’ =o +Kk, K, =sign(o,-9,), 
where K, = sign(o, —9,)|,_. - If AT >0, then «, =—1, if AT <0, then «, =+1. 

If the region c, <p <b remains elastic, then the values A, B and radius of the elastic-plastic boundary c, are 
determined from the conditions of continuity of stresses at the elastoplastic boundary, and the boundary condition 
o|,=0. So, if A and B are determined only from the conditions of continuity of stresses at the elastoplastic 
boundary, then, the following expressions take place: 


2: 
At Oe, Ba ee. (11) 
a) 3) 3¢, = 2 





The equation to calculate c, will have the form: 


3 2 
2,k) In 4). 3-26 fp pare) eer (12) 
a} 3 3b 3c, 3b Ob 


If A and B are determined from the conditions of continuity of stresses at the elastoplastic boundary, and the 








conditions o|,,=0, then the following expressions take place: 
2 3 2 
fi) es eee Ge (13) 
3b° )b ~~ 3b 3 3 
The choice of formulas (11) or (13) affects the steps of the algorithm for solving the problem, but does not affect the 
final results. 








A second plastic region will be generated at the boundary p = b if the following condition is met: 
(0, —6,) L2= kk, K, =k. (14) 
To determine the value AT = AT, when the condition (14) is satisfied, it is required to combine the system of equations 
(12), (14). Since the parameter B enters the equations (13) and (1) linearly, it is possible to obtain a separate equation 


for determining the radius of the elastoplastic boundary: 
2xykin| S.) +O aNG —P) —(0 


9 


a 3(b+¢, be, (15) 
as well as the formula for calculating the parameter Bp : 
2k, («,b° —«,c, )(b—a) 
B =p, = — (16) 


(b° —c’)ab 
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Under further loading, when the inequality B >, is satisfied, the sphere region c, <p <b goes into a plastic 


state at the boundary p=). 
A hollow sphere when exposed to heat. Consideration of material hardening. Consider the case when the 
plasticity condition has the form: 
(6,5), €&,€) )=|6, —9, —d(e; —Ef J =k. (17) 
If there are no residual deformations in the sphere before loading, then as a result of thermal heating, the plastic 
zone will be generated at the inner boundary p=a, when the conditions (15), (16) are met. With further loading, 


plastic region a<p<c, is formed. To find the stresses in this region, it is required to get the corresponding equations 


from the system of equations (3-6), (8), (17): 














2 d°o, do, 6 dabEaAT 
0 5 + 4p — —__——_| k + —————_ |= 0, 
dp dp 1+3d0-v) o(b-a) (18) 
p do, 
Oo, =— O.. 
2 dp 
The solution to the system (18) is written as: 
6, = [thin teen |S 7 
+ —v a 
‘ (19) 
_ 1 Lop _ 36abEaAT ree 
1+36(1-v) -2p(b—a) a) 


The values C,, C,, included in the formulas (19), are determined from the boundary condition o, |_.=0 and 
the condition for the absence of plastic deformations at the elastic-plastic boundary p=c;: 


C, =2«,(k-k,)c -3EaN,c, , 


2 — E : 
C, =2«,k, ing AIG 5 FON 32) (20) 


3a a a 


7 dabAT 
> +38(1-v))(b-a) 
As a result of substituting (20) in (19), we get: 


— 1 retin! 4 — ae 
* -14+36(1-v) a 


_( 2K,ke; fi dabc; EaAT \f 1 1 ri 2«, ke, _ 1 
3 b-a ap 3 la py 

Oo, = a Kk + ee ns 
1+306(1-v) 


a b-a 
Kkce’ dabc EaAT \( 2 1 Khoo ( 2 1 
- +$ ——+ ——— || —+— | |+—| —+— |. 
3 2(b-a) ap 3 la p 


From the solution to the elastic problem, it follows: if AT >0, then «x, =—1; if A7’ <0, then «, =+1. 





a 20 





Accounting for plastic compressibility. Consider the case when the plasticity function is linear with respect 
to the components of the stress tensor: 
G(20, +9, )+K(o, —9,)+N(0, —9,) =k, 1-x7), 


Kk = sign(O, —9O, ). 


(21) 


The plasticity condition (21) can be represented as: 
ao, +Bo, =kU-xT7), 
a=2¢6+K-N, 
B=¢-K+M,. 
Taking into account the introduced notation, to determine the stresses in the plastic region, we obtain the 


problem: 
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ao, +Bo, =kI- x7), 





do, 
p+ (9, -9,)=0 (22) 
6, | a= DP, 


The solution to the problem (22) has the form: 


: (M-l)k, 2Nk, (2) ~ -1 2N | 
o =|—p, +—*+—_ "+ + —__ ** ]}|- —| —_+_+——_4+_|f,, 
: ° a+p a(a+2B) /\p a+B p(a+t+2B) 


iB (M,-l)k, 2Nk, \(a\y*° (M,-1 WN, 
o, =-—| —p, +—*>*—_ + —*—__ | - — + ——*— |k,, 
(oF a+ a(a+ 2B) )\ p a+B p(a+2p) 
aAT N= abAT 
b-a’ b-a- 


To get the correct result from (20), when, for example, a+f=0, it is required to perform a limit transition 


(23) 











where the notation is introduced: M =T, — M, = YM, N=. 


when solving (22). It is easier to obtain the correct result directly in (23), while taking into account that the condition 
a+ =0 is met. In this case ¢ =0, so we get: 


2 
Go, = Ke ‘aan [4-1], |p , 
Oo a\p a 


ua d- )0-in®) +, (2-2) 9, 
a p 


as 
Ou a 





During the loading process, the plastic zone originates at the boundary p =a, when the following condition is 
met: 
AT =kt=1,)| bar SSPE ONES) 
2(1—v)(a +ab+b ) 
Research Results. Fig. 2 shows stress graphs and stress vector hodographs, when the sphere region 
corresponding to the condition a<p<c,, is in a plastic state, and the sphere region corresponding to the condition 


c, <p <b, is in an elastic state. 


Oo 


» ea Oo, =—,e€a 








M45 
a) b) 


Fig. 2. Stress graphs (a) and stress vector hodographs (5) 
for parameter values:: k = 1; v = 0.3; a=0.5; b=1; AT =170; c, = 0.57 


Fig. 3 shows stress graphs and stress vector hodographs, when the sphere regions corresponding to the 
conditions a<p<c, and c,<p<)b, are in a plastic state, and the sphere region corresponding to the condition 


c, <p <c,, 1s in an elastic state. 
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Fig. 3. Stress graphs (a) and stress vector hodographs (5) for parameter values: 
v = 0.3; a = 0.5; AT = 270; c, = 0.62; c, = 0.88 


Fig. 4 shows stress graphs and stress vector hodographs, when the sphere region corresponding to the condition 


a<p<c,,is inaplastic state, and the sphere region corresponding to the condition c, <p <5, is in an elastic state. 


fe) 
—,ea 
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a) b) 


Fig. 4. Stress graphs (a) and stress vector hodographs (4) for parameter values: 
v = 0.3; a =0.5; AT = 215; c, = 0.58 


Fig. 5 shows stress graphs and stress vector hodographs, when the sphere regions corresponding to the 


conditions a<p<c, and c,<p<)b, are in a plastic state, and the sphere region corresponding to the condition 


c, <p<c,, 1s man elastic state. 
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a) b) 
Fig. 5. Stress graphs (a) and stress vector hodographs (5) for parameter values: 


v = 0.3; a=0.5; AT = 270; c, = 0.61; c, = 0.86 
Fig. 6 shows stress graphs and stress vector hodographs when the sphere regions corresponding to the 
conditions a<p<c, and c,<p<)b, are in a plastic state, and the sphere region corresponding to the condition 
c, <p <c,, is in an elastic state. 


,ea 


. 
1 





a) b) 
Fig. 6. Stress graphs (a) and stress vector hodographs (6) for parameter values: 


v =0.3; a=0.5; AT = 79; ¢ = 0.1; n = 0.1; Ea = 0.012; y = 0.0017 


Discussion and Conclusions. The calculation results show that in this formulation of the problem, when there 
is only thermal action, the hollow sphere does not completely go into the plastic state (Fig. 2-6). Hardening causes an 
increase in the equivalent stress in the plastic region and a decrease in the radius of the elastoplastic boundary (Fig. 4, 5). 
The elastic region cannot completely disappear under loading. Plastic compressibility and the dependence of the plastic 
limit on temperature have a significant effect on the stress state of the hollow sphere (Fig. 6). 
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